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The parametric instability contribution to the largest Lyapunov exponent Ai is derived for a 
mean-field Hamiltonian model, with attractive long-range interactions. This uses a recent Rieman- 
nian approach to describe Hamiltonian chaos with a large number N of degrees of freedom. Through 
microcanonical estimates of suitable geometrical observables, the mean-field behavior of Ai is an- 
alytically computed and related to the second order phase transition undergone by the system. It 
■ predicts that chaoticity drops to zero at the critical temperature and remains vanishing above it, 

with Ai scaling as N~3 to the leading order in N. 
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I. INTRODUCTION 
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The largest Lyapunov exponent Ai is a good quantity to measure the degree of chaoticity of a generic non-integrable 
Hamiltonian system. However its numerical computation requires to compute also the microscopic dynamics for, 
sometimes, very long and, theoretically, infinite time. This may obviously turn rapidly difficult to tackle and much 
effort has been devoted to derive some asymptotic scaling laws |jj and, more recently, to get analytic estimates by 
relating microscopic dynamics with statistical averages, provided the number N of degrees of freedom is large enough 
^-^|. This latter way of analytically computing Ai as a function of e = E/N, the energy per degree of freedom, has 
proved to be remarkably efficient. It reformulates Hamiltonian dynamics in the language of Riemannian geometry, 
q ■ using the fact that the natural motions can be viewed as geodesies of a suitable Riemannian manifold Q. Chaotic 
motion then reflects into the instability of the geodesic flow which depends on curvature properties of the manifold. 
• This geometric formulation of the dynamics has been known for long and led to fundamental results in abstract 
' ergodic theory when the ergodicity of geodesic flows on compact manifolds of negative curvature was demonstrated 
by Hedlund and Hopf in 1939, and later exploited by Krylov ||. However, when more physical Hamiltonian systems 
come into play, such as coupled nonlinear oscillators, a major source of chaos appears to be parametric instability 
activated by a fluctuating curvature along the geodesies, even when curvature is always positive This has been 

exploited in the theoretical model proposed by M. Pettini and coworkers. Modeling the effective curvature felt by a 
geodesic by a gaussian stochastic process, with mean the average Ricci curvature and variance its fluctuations, and 
under the ergodic hypothesis replacing the previous geometrical quantities with their averages kq and a\ according 
to the natural ergodic measure, i.e. in the microcanonical ensemble, they derive the following expression for Ai 
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and r, a timescale for the stochastic process estimated as: 



(3) 
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In this article, we apply these geometrical tools to a mean-field Hamiltonian system of globally coupled rotators 
exhibiting a second order phase transition at a certain critical energy e c . We analytically estimate the parametric 
instability contribution to Ai(e) and predict a neat distinction between the two cases: e < e c and e > e c . Numerical 
simulations |^|,|o| seem to qualitatively support the analytical conclusions. The remarkable behavior of the Lyapunov 
exponent in the mean-field limit, as a consequence of the simple expressions of relevant geometrical quantities as 
functions of the order parameter, could then be a dynamical signature of the phase transition. 

The model at hand will be described in Sec. II and some useful geometric expressions derived there. A detailed 
derivation of the largest Lyapunov exponent Ai as a function of the energy density e will be exposed in Sec. Ill, Sec. 
IV being devoted to comments and conclusions. 



II. THE MEAN-FIELD MODEL AND FIRST USEFUL GEOMETRIC EXPRESSIONS 



Here we study the so-called mean-field Hamiltonian X-Y model, which can be considered as a toy-model for inves- 
tigating long-range interactions in Coulomb systems [ pd| , ^2| . The dynamics of N interacting particles moving on the 
unit circle II = [0; 2tt] derives from the following Hamiltonian 
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[l-cos( qi -q r )] = K + V(q) 
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where K and V stand for the kinetic and the potential energy respectively. Constant c may be rescaled to +1, or —1 
by a change of variables. The scaling factor for the potential energy ensures that the interaction energy is extensive 
and emphasizes its mean-field nature. Thus, in the following, we would not deal with the usual thermodynamic limit 
with fixed density, but rather with the mean-field limit N — > oo, ■§> — > e, e finite. Note that the total momentum is 
also a constant of the motion. However this will not affect following calculation since the potential only depends on 
positions. 

The equilibrium statistical mechanics of this model can be exactly derived [ft3| . In the case of an attractive potential 
(i.e. c > 0), which will be assumed in the following, that is in the ferromagnetic-like case, it predicts a second-order 
phase transition with order parameter ||M|| where M is the mean-field magnetization-like variable defined as: 
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This phase transition can be easily conjectured by observing that at small energy ||M|| = O (1) with a clustered phase, 
whereas at large energy, the central limit theorem predicts that ||M|| = 0(N~i) with particles having random ballistic 
motions. It is also interesting to note that introducing the global variable M enables to re-express the equation of 
motion of any particle as: 

iji = — c ||M|| sin(<7i — <p) where 4> = arg(M) (6) 

that is the equation of a perturbed pendulum, -the full system being closed by adding the evolution equations for 
||M|| and (j). 

Let us now first express in the framework of the Eisenhart metric the Ricci curvature associated to this system, 
then derive the microcanonical averages of the geometrical quantities involved, via the canonical ensemble which leads 
to simpler calculations. Recall here that in the limit of infinite size, that is N — > oo, the averages of thermodynamic 
observables in different ensembles coincide H|, but not their fluctuations flEf] . Therefore, in order to get the fluctu- 
ations of an observable / in the microcanonical ensemble, it will be necessary to add a corrective term according to 
the formula derived in |16| , which is not valid at the critical point: 
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So, with the Eisenhart metric, Ricci curvature reads K R (q) = AV, where A stands for the euclidian Laplace 
operator in the configuration space, so that the average Ricci curvature ||, defined as k R (q) = , is 

Moreover, a straightforward calculation gives: 
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n<z) = ^(l-HM|| 2 ) (10) 

Thus we obtain the key expression that the mean Ricci curvature reads simply in terms of the order parameter, the 
mean-field magnetization M as: 

k R = c\\M\\ 2 (11) 

up to a O (iV -1 ) term, which, as far as the mean-field limit is concerned, gives a vanishing contribution and will be 
ignored. It will only play a part in corrections above the transition. It should be pointed out that this expression 
for the mean Ricci curvature as a smooth function of the natural order parameter, the magnetization, is not claimed 
here (since not proved) to be a generic property of, for instance, some class of mean-field Hamiltonian systems. At 
present we should thus consider the results obtained in this article as peculiar features of the model at hand. As 
only positions-involving quantities come into play, let us now focus on the contribution of the potential energy to the 
partition function in the canonical ensemble at temperature T = (3~ 1 (with k R = 1): 

ZM = V *M-PV(q))d N q = exp(-/3^) /„„ eMPf \M\ 2 )d N q 

Then, using the integral representation of gaussian functions, we get: 
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Z c ((3) = exp(-/3^) J nN I [J R2 exp(-u 2 + 2^(3^u.M)du 

= exp (-/3<f ^ r K2 exp(-u 2 ) [I (2^\\u\\)] N du 
= (2Tr) N f-J™rdrexp(-Ni;(r,P)) 

where ip{r, f3) = — ln(io(r)) + 4p and where /„ stands for the modified Bessel function of order n. 
Then, according to the saddle-point method, in the limit N — > oo the previous integral is fully dominated by the 
minimum of ij) obtained by solving the consistency equation d r ip(r,[3) = that is: 

pc I {r) 

When /3c < 2, tp is minimal for r = 0, which corresponds to a vanishing magnetization. For /3c > 2, ( |l2| ) admits a 
non- vanishing solution noted r*(/3), the phase transition taking place for [3c = 2 i.e. for T c = c/2 and s c = 3c/4. 
Before examining these two cases, we establish some useful canonical relations: as (V(q)) c = —dpln(Z c ) and 

(V(q) — (V^q))) 2 ^ = 9 2 ln(Z c ), one obtains respectively: 

(k B ) c = c+^d P HZ c ) (13) 



(6*K R ) c = i {{K R {K R )f) c = la 2 \n(Z c ) (14) 
Moreover the energy density e{(3) is given by: 

e(/3) = h~^ dp{lnZc) (15) 

In the following, when dealing with microcanonical estimates, this expression will be implicitly systematically used 
to express f3 as a function of the energy density. We define also the two notations = (k R )^ and cr 2 = (5 2 Kr/ . 
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III. ANALYTIC ESTIMATE FOR Ai BELOW AND ABOVE THE TRANSITION 



Let us now derive the analytic estimate for Ai below and above the transition. Below the critical energy, the 
saddle-point method gives: 



As the ensemble averages (k R ) c and (k R ) coincide in the mean-field limit, this gives: 
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That is: 



r*(/3) 2 



Remember that is proportional to the square norm of the magnetization so that we expect it to exhibit the 
same behavior at the transition point with twice the characteristic exponent. Actually, a straightforward expansion 
near the transition leads to 

2Q3c-2) 8 
jg = YTI~c i£c ~ S) £c ~ £ - 

Taking into account the correction (Q) and noting that dp (k R ) c = \ (b^K^ , one finally obtains: 

(6*K R ) =^K R ) (l + t^ KR ) V 1 (is) 



with <<5 2 A^) c - ^ln^-i^r* - £). 

Figure 1 displays the behaviors of both the average Ricci curvature kq and fluctuations a Kl with the control 
parameter c put equal to 1 in both figures. Using (|l7],[l8|), one can then derive Ai(e) in the clustered phase. The 
result, obtained through (0J|J|), is reported in Fig. 2. When e approaches e c , expanding the expression for the largest 
Lyapunov exponent Ai(e) provides the scaling law 

Ai(e)cx(e -e)* (19) 

associating thereby a critical exponent, equal to 1/6, to the dynamical observable Ai. 
Above the critical energy, one obtains in the same way: 

Z c (P)^(2n) N eM~N^)(l-^j (20) 

Here, as ||M|| 2 becomes of order O (iV -1 ), we shall use the full expression k R = c 

ll M H 2 - § +0 (N- 2 ). Then: 

i.e. the microcanonical average of the Ricci curvature vanishes in the mean-field limit. 
Similarly: 

(5 2 K R ) c = ^dpn(Z c ) = ^ (2- (3cy 2 = O (N^) 

As e(/3) ~ m + § , the correcting term needed to get the microcanonical fluctuations is of order TV -2 , thus negligible. 
And then: 

(5 2 K R )^ ~ ^ (2 - (3c)- 2 = O (TV- 1 ) (22) 

We can keep in further calculations the dominant order in N, and derive the scaling law with TV for the largest 
Lyapunov exponent. Using expressions (^J^J^) , in the limit N — > oo, one obtains 

(2 - (3c) 2 
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IV. COMMENTS AND CONCLUSIONS 



Let us first comment here on the reliability expected for the expressions just derived. As developped in refs. ^^), 
the geometrical approach aims at extracting information on, at least, an average degree of chaoticity of the dynamics 
from mean global geometrical properties of the Riemannian manifold constructed from a given Hamiltonian. This 
implies the crucial assumption of ergodicity, as a way of bypassing the knowledge of the trajectories, i.e. the numerical 
integration of the equations of motion. This ergodic hypothesis is not expected to be realized in the integrable limits 
of small and large energy, the latter following from the boundedness of the potential energy in (^). However, it is well 
known that chaos is not a necessary condition for ergodicity, the most striking piece of evidence being provided by the 
ideal gas of point particles, for which there is no velocity mixing at all. Also recent studies jlq] have emphasized that 
crgodic-like properties should depend mainly on the observable at hand, irrespectively of the degree of chaoticity of the 
dynamics. Concerning our model, S. Ruffo already observed in jl^] a good agreement between Gibbs predictions and 
numerical simulations for the observable ||M||. Moreover, in the mean- field limit, this happens even in the integrable 
limit of large energy, an explanation for this being provided by a result of M. Kac |T^,[l9|], so that the mean- field 
magnetization appears like a good observable with respect to ergodicity. Therefore it is not surprising to observe that 
numerical calculations of the mean Ricci curvature and its variance fit well the microcanonical predictions presented in 
Fig. 1 ||,^(| , except in the vicinity of the phase transition where finite- AT effects dominate. Concerning the transition 
region, as noted before, the formula pl| used to get fluctuations in the microcanonical ensemble from canonical ones 
is not valid at the critical energy. Therefore we should exclude in our conclusions a small neighborhood of e c , all the 
smaller as N is large. So the analytic estimate for Xi(e) in the mean-field limit is expected to be quite reliable except 
maybe for small e and in the vicinity of e c . It should also be noted that the timescale r estimated as (|J), that is the 
time under which the effective curvature felt by a geodesic cannot be regarded as a random process, is the less solid 
point of the geometrical modeling as (||) relies mainly on phenomenological arguments. Then it can, if necessary, 
be slightly adjusted to fit numerical calculations. Nonetheless that estimate for r is also a powerful tool, as it provides 
a natural timescale, depending on e, that should be taken into account to connect for instance results for mappings 
|jj to results for continuous flows as it is the case here. 

Keeping these remarks in mind, we can now comment on the results obtained in Sec. III. Expression ( p3| ) means 
that, in this mean-field model, above the critical energy, chaos does not survive to the limit TV — > oo. This can 
be conjectured straightforwardly from the equation (^|) governing the time evolution of any particle, which predicts 
ballistic motion as ||M|| vanishes above e c . Moreover, one obtains the scaling law N~3 for the largest Lyapunov 
exponent to the leading order in N. The same scaling law has been found numerically by Latora et al. 0|. A rather 
nice fit (see Fig. 2) is also obtained with Yamaguchi's simulations jlO| on a wide range of e, except in the vicinity 
of e c , where finite size effects smooth the transition. Here strong metastability related to critical slowing down may 
also affect numerical results with relaxation times towards equilibrium increasing greatly with N. Besides, for a given 
N large enough, expression ( p3| ) rightly gives a vanishing Lyapunov exponent in the integrable limit of large energy 
where rotators tend to behave as free particles. 

Concerning the transition region, in spite of the above mentioned remarks on the validity of our results at the critical 
energy, let us mention the remarkable features exhibited by Figs. 1,2: kq, <t k and Ai display singular behaviors at the 
critical point. Here curvature fluctuations exhibit a discontinuity which is similar to the "cusp" numerically observed 
in j|. In our case, this appears as a direct consequence of the second order phase transition exhibited by the model 
and, following equation (|ll|), of the expressions of the different parameters used in the geometrical approach in terms 
of smooth functions of the order parameter. Following conjectures exposed in ||, the geometrical meaning of these 
singular behaviors might be that a topology change of the "mechanical" manifold underlying the dynamics occurs at 
the critical energy. 

Finally, as for Ai, its maximal value would be reached slightly below the critical point and not at the critical point. 
Numerical simulations made in || for 20000 particles show such a tendency. Moreover, when e approaches the critical 
energy, calculations ( |l9| ) show that Ai goes to as (e c — e)s . This suggests that a critical exponent could be associated 
to the largest Lyapunov exponent as a dynamical observable. 

Further studies should inspect more precisely the region where the amplitudes of the curvature and fluctuations are 
comparable, around e = 0.45 (see Fig. 1). As observed in other models, for such a situation strong stochasticity may 
be expected. A more refined treatment may imply some corrections to the gaussianity of the effective curvature, that 
would take into account further moments of the mean Ricci curvature. Also the vicinity of the critical energy, as 
well as a possible extension of the results obtained in this article to a larger class of mean-field Hamiltonian systems 
deserve obviously further investigations. 
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Figure captions 



Fig. 1: Analytic expressions for the microcanonical averages of the average Ricci curvature, kq (solid curve) and 
of its fluctuations a K (dot-dashed curve) in the mean-field limit, below and above the phase transition. 



Fig. 2: Analytic expression for the largest Lyapunov exponent Ai in the mean-field limit (solid curve) below and 
above the phase transition. Analytic corrections (dot-dashed curves) to mean-field limit for finite N with N = 80 and 
N = 200 above e c . Here the derivation does not restrict to the leading ter m ^23| ) but computes (d,||,§) up to further 
orders, as N is not very large. There is a nice fit with results exposed in [|l 0f| apart from the vicinity of the critical 
energy. 
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